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ABSTRACT 



The identity of dark matter is one of the key outstanding problems in both particle 
and astrophysics. In this thesis, I review some candidates of dark matter, especially 
WIMPs (weakly interacting massive particles) which is one of the best candidate so 
it is called that WIMPs miracle. In addition of this, there are also some theories of 
modification of gravity, by changing the law of gravity, it could be possible that the 
dark matter observations are explained. Until the dark matter particle is detected, there 
is some room for uncertainty. So we should consider every part of the problem and 
solve it. Dark matter problem is covering a large area so every possibility is important. 
So f(R) gravity is also reviewed in this thesis and some theories are considered as a 
possible solution of dark matter problem. Finally we highlight that, even in the case of 
WIMPs or another particles solution, f(R) gravity is also can be used for this problem. 
However, last words will be said by experiments. 



Keywords: Dark Matter, Weakly Interacting Particles(WIMPs), f(R) Gravity, Cos- 
mology. 
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Karanlik Maddenin gizemi hem astrofizik hem de parcacik fizigi acisindan 50k onemli 
bir problemdir. Bu tezde, karanlik madde adaylanni, ozellikle WIMPs mucizesi olarak 
bilinen ve en giiclii aday olan WIMPs (zay if etkile§en kiitleli parcacik)i inceledim. 
Karanlik Maddenin gizemini cozmeye caligan bundan ba§ka teoriler de vardir, ornegin 
gravitasyonun modifiye edilmesi gibi. Karanlik madde parcaciklan deneylerde bulu- 
nana kadar, bu gizem siirecektir. Biz bu gizemi her acidan inceleyip coziim iiretmemiz 
gerekir ciinkii cok geni§ bir alani kapsayan bu alan icin her bir ihtimal bile cok onemlidir. 
Karanlik maddeyi a§iklamaya aday olan f(R) gravitasyon teorisini de bu tezde in- 
celedik ve coziim iiretmeye §ali§tik. Sonu5 olarak, gerek WIMPs olsun gerekse baska 
par5aciklar olsun, f(R) gravistasyon teorisi de giizel bir alternatiftir ve kullanilabilir. 
Tabii hangisinin dogru oldugunu deneyler onaylaycak. 



Anahtar Kelimeler: Karanlik Madde, Zayif Etkilegen Kiitleli Parcaciklar(WIMPs), 
f(R) Gravistasyon Teorisi, Kozmoloji. 
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Chapter 1 



INTRODUCTION 

1.1 Observational Evidences for Dark Matter 

Dark Matter is a very strange form of matter. It does not interact with known objects 
because of the lack of electromagnetic interaction. However, it has a mass and we 
know that it exists. The main reason for this view is from gravitational interactions on 
the visible matter which holds galaxies and galaxies of clusters together. The idea of 
the existence of a new form of matter escaping electromagnetic detections, traces back 
at least to 1915. The name "dark matter" was firstly used by the Jacobus Kapteyn in 
1922 in his studies of the motions of stars in our galaxy [2J. He found that there is no 
need of Dark Matter within and around the Solar System. In 1923, Jan Oort arrived to 
an opposite conclusion that there should be twice as much dark matter as visible matter 
in the Solar vicinity. This is the first claim of evidence for dark matter. However, since 
later observations have disproved this early claim, in 1933 the discovery of dark matter 
is usually credited to Fritz Zwicky who made the first correct claim about the existence 
of dark matter.He concluded that the velocities of the galaxies in the Coma cluster 
are much larger than the expected one. Zwicky found that orbital velocities are almost 
a factor of ten larger than expected from the total mass of the cluster [fT6l . Zwicky 
concluded that, in order to hold galaxies together in the cluster, the cluster should 
contain huge amounts of some Dark (invisible) matter. The galaxy rotation curve v, 
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where r marks the distance from the galaxy's center. 



^ = <££, (i.D 

r r z 



where M(r) is the mass inside radius r and = 6.7 x 10 8 cm 3 g [ sec 2 . After rear- 
ranged the above equation reads 




(1.2) 



If the mass is only inside the radius, the velocity should decrease when there is 
large value of radius as 

vo-l/v^ (1.3) 

, where G = 1 . 

However, for values of r bigger than the size of the galactic bulge, the velocity 
of stars remains approximately constant! 15] . It means that there is a some unknown 
energy there such that 

M(r)oc r (1.4) 
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Figure 1.1. 

Distribution of dark matter in NGC 6503. Rotation curve of NGC 6503. The curves 
labeled gas, disk, and halo represent the contributions from the mass in gas, stars, and 
a dark halo. This plot of the rotation curve is from Bertone et al. (2005b) using data 

from Begeman et al. (1991) EM 
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For an orbit around a compact central mass, such as planets in the Solar system, we 
get Kepler's third law which is given at Eqn.(1.3). 

As an example if the star is placed in the center of a galaxy, that is the mass density of 
a galaxy decreases as a power law i.e 

poc r - n (1.5) 

The mass inside radius r is 

M(r)oc I r 2 r - n drocr 3 - n (1.6) 



Thus the rotation velocity is 



v(r)=r 1 -"/ 2 (1.7) 



Near the center of the galaxy, observed rotation curves of a galaxy v(r) increase with r 
from Eqn.(1.3), however then suddenly becomes v(r) constant. From the Eqn.(1.5), 
the density profile is 



pocr~ 2 (1.8) 



On the other side, the density of stars falls faster when goes to the edges of the galaxy. 
The total mass from stars, gas and dust clouds, is very very small so that they are 
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not able to provide the rotational speed at large distances. This means that there is 
another invisible, massive matter. This matter would be only dominant in the outer 
parts. The dark component forms a dark halo which holds galaxies and galaxies of 
clusters together. 

At galactic scales, there is also additional evidence for dark matter which is based on 
the mass modeling of the detailed rotation curves [|4]|. An evidence of Dark Matter 
from cluster of galaxies was first noticed, as mentioned already, by in 1933. On the 
cluster scales, the Q-m is today's mass density including baryonic mass and dark matter 



As known that matter relic density component is bounded by the observation of the 
large scale structure. As reported by Sloan Digital Sky Survey (SDSS) Q.Mh 2 = 0.135 
for £l b h 2 = 0.025 0. The mass density of baryons and the mass density of matter in 
the universe by using WMAP data are found to be |[T3~ll8l 



ttHil 



£2 M ~ 0.2 -0.3 



(1.9) 



£l b h 2 = 0.0223 



(1.10) 
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Figure 1.2. 

WMAP data. Credit: NASA / WMAP Science Team £D 
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Figure 1.3. 

CMB Temperature fluctuations. Image from NASA / WMAP Science Team flU. 
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Figure 1.4. 

WMAP observes the first light of the universe. Credit: NASA / WMAP Science Team 

m 

a M h 2 = o.m (l.ii) 



From WMAP's data and large scale structe observations 's data for relic abundance 
are consistent. In physical cosmology, Big Bang nucleosynthesis provides us also a 
constraint lfTOl 

The difference between Q. M h 2 and £l b h 2 is the clue of the existence of the non-baryonic 
dark matter. 

There are many experiments which support the existence of Dark Matter such as CMB 
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anisotropy, galaxy surveys, etc. so we think that there is Dark Matter. Unfortunately, 
the nature of dark matter is now an open problem. So this is a very hot topic in physics. 

Furthermore, there are some additional constraints on the amount of dark matter in the 
Universe. We know very well from the concordance cosmology measurements(it's ba- 
sically a combination of looking at the microwave background fluctuations, the baryon 
acoustic oscillations, and the supernova brightness versus redshift relation)what are the 
relative fractions of dark energy and mass in the Universe. We also have pretty good 
constraints from primordial abundances of the elements on what the relative fractions 
are for baryonic and non-baryonic dark matter because if there were too high a density 
of protons in the early universe, then there would have been more fusion in the early 
universe irrTll . 

1.2 Baryonic Dark Matter(BDM) 

MACHO(massive compact halo object); planet-like objects, plasma and other com- 
pact objects are the candidates for BDM. Other examples of this class of dark matter 
candidates include primordial black holes created during the big bang, neutron stars, 
white dwarf stars and various exotic stable configurations of quantum fields, such as 
non-topological solitons. Gravitational Microlensing is one of the ways to detect these 
objects. However, this method only gives us the information about the mass, distance 
and velocity of them and nothing more. 

Heavier, relatively faint objects such as old white dwarfs and mini black holes,are 
the most important candidates for this BDM. If there are primordial mini black holes 
produced in the big bang before BBN(Big Bang Nucleosynthesis), they would not con- 



9 



tribute to £V 

0.07A/ Q is required by a star to start to shine as a star by igniting thermonuclear fu- 
sion. Brown dwarfs are not completely dark because they radiate faint thermal radia- 
tion. However, the typical mass of these MACHOs should be 0.5M q , which is much 
larger than the brown dwarf mass. In addition, white dwarf has a mass which is suit- 
able for dark matter but there is not enough number of observed white dwarfs. Hence, 
BDM is not as good candidates for dark matter as nonbaryonic dark matter. 
1.3 Nonbaryonic Dark Matter 

Hot Dark Matter (HDM) and Cold Dark Matter (CDM) is considered as a nonbaryonic 
dark matter. Other more special options have been considered, such as for example 
warm dark matter (WDM). 

Dark matter particles that are called HDM, are relativistic at the time of matter-radiation 
equality time. This requirement implies that their mass has to be lower than about 100 
eV. Their thermal velocities, when structure formation begins are large because of the 
small mass. In contrast, CDM has large mass so its velocities are negligible. Candi- 
dates between hot and cold are called warm dark matter. 

Effect of the pull of gravity causes to form galaxies and galaxy clusters from early 
homogeneously distributed. This is called structure formation. Today the best way 
to differentiate between HDM, WDM and CDM is through the observed large-scale 
structure, i.e. CMB which shows the seeds of structure. We show in Fig. (1.3) the 
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Figure 1.5. 

Comparison of the expected halo of the Milky Way and the galaxies M3 1 and M33 in 



CDM and WDM models. From http://www.clues-project.org/images/darkmatter.html 



results of a simulation of the halo of dark matter around the Milky Way and two other 
galaxies. For CDM, there is more substructure and satellites around the galaxy, while 
their formation is suppressed for WDM. According to observations, there is an order 
of magnitude less satellites observed around the Milky Way than predicted by CDM 
models. However, the observations are not complete, and the discrepancy may also be 
due to other causes than WDM. 

The most common candidates for non-baryonic dark matter are Weakly Interacting 
Massive Particles, or WIMPs flU. Cold dark matter has a most important candidate 
which is WIMPs. Because they have neutral charge and they are long lived parti- 
cles. Only weak force and gravitational force effect them. Large amount of WIMPs 
are produced and then the interaction rate of WIMPs was decreasing so decoupling 
increasing between them in the thermal bath of the early universe when cools down, 
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like neutrinos, but are much heavier, so that they are a form of CDM. After the de- 
coupling, their relic density remains constant as consistent with experiments. This is 
called WEVlPs miracle. Furthermore the interactions of some dark matter candidates 
are stronger or weaker. For example, gravitinos have only gravitational interactions, 
while TIMPs (Technicolour Interacting Massive Particles) have much stronger interac- 
tions. Neutralino (bino+wino+2higgsino) lightest supersymmetric particle is the best 
candidate for WIMP in Minimal Supersymmetric Standard Model (MSSM)|HH E3- 
The other candidate for WIMP is the lightest Kaluza Klein (LKP) which is coming 
from Universal extra-dimensions (UED) model . 
1.3.1 Hot Dark Matter (HDM) 

For the HDM neutrinos are candidates. Energy density today is dominated by their rest 
masses so to make a significant contribution to the density parameter, neutrinos should 
have a rest mass about 1 eV. The present upper limit to HDM in the form of massive 
neutrinos leads to the conservative mass limit from WMAP data. 



£m v <0.62eV (1.12) 

Therefore the maximum contribution of neutrinos is Q. v h 2 < 0.02, at most the same 
order of magnitude as baryonic matter. If neutrinos were the dominant form of matter, 
there would be a lower limit on their mass from constraints on the phase space density, 
called the Tremaine-Gunn limit. Essentially, in order to achieve a certain rotation 
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velocity for galaxies, you need a certain amount of mass inside a given volume, and 
the Pauli exclusion principle constrains the number of particles you can pack inside a 
given volume. Even though we know that neutrinos are a subdominant component of 
dark matter, the Tremaine-Gunn limit applies to any fermionic dark matter candidate, 
even if its distribution is not thermal. (There is no similar lower limit on the mass of a 
bosonic dark matter particle.) 
1.3.2 Cold Dark Matter (CDM) 

Most of the matter in the universe has to be CDM because of data on large scale struc- 
ture combined with structure formation theory. WIMPs (Weakly Interacting Massive 
Particles) are the main candidates of CDM. Dark matter candidates are with masses 
in the 100 GeV to TeV range. Particles with GeV to TeV masses and weak-scale in- 
teraction cross sections are known as Weakly Interacting Massive Particles (WIMPs). 
There are also another candidates for DM. However, WIMP dark matter is compelling 
for two main reasons: 

1 . In the early universe, when WIMPs were in thermal equilibrium annihilated with 
each other with a weak scale cross section. As the universe expanded and has 
frozen out, WIMPs cooled to obtain a relic density today which is equal to the 
measured dark matter density. 

As mentioned before that this is known as the WIMP Miracle. 

2. The other reason is that WIMP dark matter are present in many of the favored 
extensions of the Standard Model. For example, the Minimal Super symmetric 
Standard Model (MSSM) has a favourite WIMP candidate. 
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Chapter 2 



EXTENDED GRAVITY: f(R) GRAVITY 
2.1 Introduction to f(R) Gravity 

By changing the law of gravity, it could be possible that the observations are explained. 
Until the dark matter particle is detected, there is some room for uncertainty. The 
problem for modified gravity is that there are so many different observations for dark 
matter, in different physical systems: motions of stars in galaxies, motions of galaxies 
in clusters, gravitational lensing, large-scale structure, CMB anisotropics and so on. 
Gravity has to be adjusted in a different manner for these different observations, and 
the resulting models are somewhat contrived. Expressed another way, the dark matter 
scenario is very predictive: the simple hypothesis of a massive particle with weak 
couplings to itself and to the Standard Model particles explains a number of disparate 
observations and has made successful predictions. 

In 1919 and 1922 respectively by Weyl and Eddington [fT9l l20l gravitational action 
firstly studied. However, until 1960s, there was not any indications of complicating 
the gravitational action. Because of cosmological problems such as the dark energy 
problem, dark matter problem, this type of modification is needed to fix the problems. 
At f(R) theory, the main idea is to take action as a function instead of a constant 
curvature such as: 
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(2.1) 



By using the f(R) function , it provides usage of many different properties, which are 
very important for explaining unknown matters such as dark matter and dark energy. 
For instance, considering f as a series expansion i.e Il22l [231 1241 |25l |26l l2~7Tl 



/(K) = ... + | + f-2A + K + g + g (2.2, 



where a, and P ; coefficients have the appropriate dimensions. So it can be seen that a 
number of interesting terms appear. ||2TI 

Without matter action is given above, now after adding matter on the action, our f(R) 
action is 



S = J d 4 x^[f(R) +S M (g f »,V)}. (2.3) 



Varying the action with respect to provides 
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(2.4) 



where f R =% , □/* = 7^ (v^ (g)) and V^f R = g^ [(f R )^ a -^ a (f R ) 



So 



1 2 /—nfdf 



(2.5) 



0v7 d / 1 JXW f J A 1 



(2.6) 



v'v lA = v'v 1 ^=,»v 1 v 1 ^=,"( 



ll r l 



(2.7) 



(2.8) 



in which a prime means d/dr and its variation with respect to g MV yields, after some 
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manipulations and modulo surface terms, the field equation 



f R R»v ~ -^Suv = V„V V /^ - glN Uf R + K^v (2.9) 

where is the covariant derivative associated with Levi-Civita connection of the met- 
ric, and □ = W%. 

These equations above are fourth order partial differential equations in the metric. The 
alternative name "fourth order gravity" used for this theories because in the first two 
terms on the right hand side fourth order derivatives appear. 

Furthermore, the contraction of the above equation yield the trace equation 

f R R-2f + 3Df R = KT (2.10) 

where T = = 

the trace equation can be used to simplify the field equations and then keep it as a 
constraint equation. 

Therefore 

f R R^ - ^ (f R R - Df R - KT) - V v V v f R = K7? (2. 1 1) 
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The trace of this equation is 



f R R-2f(R)+3Df R = SnGT (2.12) 

where T = g^ v T^ v 

Here we can see that when T = , not anymore says that R = or constant. 
Another important thing is that when R = constant and 7^ v = 0, it reduces to 

f R R-2f(R) = (2.13) 

Maximally symmetric solutions is affected by f function, and they will be Minkowski, 
de Sitter or anti-de Sitter. Hence, We can now write the field equations in the form of 
Einstein equations as 




(2.14) 



After some algebra the above equation, takes the form 
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F(R)R, V -V^F(R) -kT, y = -g^( l -f{R) - V a V a F(i?)) (2.15) 



where F(R) = df/dR = f R . 

For solving above equation we need a metric and its components ETTl . Here we are 
using spherically symmetric metric with radial components B(r) and X(r) : 



X(r) 

ds 2 = -B(r)dt 2 + jj-tdr 2 + r 2 dQ. 2 (2.16) 



Ricci scalar of this metric can be written as 



. . 2 X' ,B' 2R 1 . „ AB' 2R 

' ?w = ^ + x5 ( 2 + y ) -x (s + - + ^ (2 - 17) 



For empty space if we solve the Einstein Equation with the form of f(R) gravity, two 
independent field equations are found : 



X' 2rF" 
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J7 1 Y 1 9 P 1 Y 1 9 

B " + (L. _ _ ± (L. _ ±-\ B - 4x = (2.19) 

V F 2X ; r y F 2X } r 2 v 7 



where the prime means d/dr. After that are choose our ansatz with constants Fq, Xq and 
Bq as 



F(r) =F r n ,X(r) = X r a ,B(r) = B r a (2.20) 



Constans are bounded by two constraints as: 



n + Z Z 4 



After that 



is obtained. 



where 
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This will be the well known solution of Schwarzschild when a or n goes to zero. After 
using those constraints and definitions, an asymptotic form of f(R) where Ao= constant 
of integration is found as 



f(R) 



/o(2*o)(5> 



(2.24) 



In the paper of A Model of f(R) Gravity as an Alternative for Dark Matter in Spiral 
Galaxies they took Fo = 1 when n is going to zero and action goes to FqR + Aq. 

For the application of this in the rotation curve, geodesic equation is used in the weak 
field approximation. The test particle is rotating around the central mass and obtains 
as 



where the dot means d/dt. After using the affine connection which we found by using 
metric : 








(2.25) 
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Br * _ £rsm(e) 2 e _ «n(9) 
1 ee - y^-^W* 1 # ^ > A # ^ — ( z - z/ ) 



velocity for a test particle is obtained as 



2 rc^i? c 2 a^ r a 
v (2.28) 



After considering some assumption for the large distances from center of galaxy, such 
as using r a = 1 + aln(r) where a is smaller than 1, one finds that 



o c aBft 

vl = (2.29) 



For taking Bq = 2^/jM ga i axy where fi is the proportional coefficient with mass inverse 
dimension and M is mass of galaxy, hence asymptotic velocity is found 
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vL = C 2 a^/jM ga i axy 



(2.30) 



where a and /j depend on experimental datas. For instance, if this values taken as 
aV = (voo/c) 4 /M galaxy = I0~ 54 kg~ l , the velocity is found exactly lOOkms^ 1 and this 
value is consistent with Tully-Fisher's relation. [12111 

Furthermore, another good testing is made by using the Modified Newtonian dynamics 
(MOND). 

0C 2 a/ = 0.98x KT 54 ^ 1 (2.31) 
for the velocity of typical galaxy. 

There are many difficulties of this theory. One example which gained a lot of publicity 
a few years ago is the Bullet cluster. A solution of the Dark matter conundrum in 
terms of a new elemantary particle is quite successful to explain the Bullet cluster, 
while explanations in terms of modified gravity models encounter serious difficulties. 
Another problem in this way of solution is that ansatz is choosen by hand and not 
naturally comes out. 
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Chapter 3 



DARK MATTER PRODUCTION 
3.1 Thermal Production 

WIMPs and many other dark matter candidates are in thermal equilibrium in the early 
universe and decouple once their interactions become too weak to keep them in equi- 
librium. Those particles are called thermal relics, as their density today is determined 
by the thermal equilibrium of the early universe. When the annihilation processes of 
WIMPs into SM particles and vice versa happend at equal rates, there is a equilibrium 
abundance. When the Universe cooled down and the rate of expansion of the universe 
H exceeds the annihilation rate, the WIMPs effectively decouple from the remaining 
SM particles. So the equilibrium abundance drops exponentially until the freeze-out 
point.The abundance of cosmological relics remained almost constant until today lTT8l . 

If the candidate is stable or has a long lifetime, the number of particles is conserved 
after it decouples, so the number density falls like a~ 3 . Specifically, we use the Lee- 
Weinberg Equation [29]. It describes annihilation and creation of % particles [|28]|30l . 

+ = " < OV > - n leqn) (3-D 
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with n^ eqn ,the number density of the relic particles in equilibrium. Furthermore Hub- 
ble constant make expanding universe effect. The thermal average of the annihilation 
cross section a multiplied with the relative velocity v of the two annihilating % parti- 
cles. The first (second) term on the right-hand side (RHS) of the equation describes the 
decrease (increase) of the number density due to annihilation into (production from) 
lighter particles. Kinetic equilibrium is satisfied by the Lee- Weinberg equation for % 

Now we use Hubble constant effect which provide universe expansion by considering 
the evolution of the number of particles in a portion of comoving volume given by 



N = n x R 3 (3.2) 



We can then introduce the convenient quantity 



- n y R 3 

N = — — (3 3) 

N(0) K } 



such that iV(0) = 1. 

In addition, since the interaction term usually depend explicitly upon temperature 
rather than time, we introduce the x = ™, the scaled inverse temperature. 

During the radiation dominated period of the universe, thermal production of WIMPs 
takes also place . For this time, one can find the expansion rate 
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H=J™g* (3.4) 



and then 



„ {x=l)=m \l^ (3.5) 



where H(x) = tf ^ 2 ^ and where x refer the number of intrinsic degrees of freedom 
such a like that spin or color for % particle. 

Hence we can write the last format of our equation a 



dN < av > N(0) 2 2 



where 



iV(0) 3^[3]m 3 



R 3 4Tl 2 x 3 



After all little tricks , our Lee -Weinberg equation can be recast as 



(3.7) 



f = -^-Nl qn) (3.8) 
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where we defined 



W(x) = §t/^^/ < <3v >= 0.055-^mM Pl o (3.9) 
4n V 4IP ^gj ^jft 



where C, is EulerRiemann zeta function To integrate the Lee -Weinberg equation, 
we need to have an expression for the equilibrium number density in comoving vol- 
ume. Once the particle is non-relativistic, the difference in statics is not important. 

The general equation of equilibrium number density in comoving volume is 



N ^Un=jnf y ^7f=2 d y < 3 - 10) 

^P\Jo l+expv* 2 +r 



For the nonrelativistic case at low temperatures T 3> m x one obtains, 



^W e ^ = ^j^ 3/2 exp- (3.11) 



At great value of temperatures, dark matter quickly is destroyed by its own antiparti- 
cle. Shortly after that, until T x = N(x) < ov H, T has dropped too much which is 
under the value of m x (T « m x ), also the number density of dark matter drops expo- 
nentially. The temperature at which the particle decouples (The time when the number 
of particles reaches this constant value) is called freeze-out temperature 7>. Hence the 
number density per comoving volume becomes almost constant because dark matter 
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particles are no longer able to annihilate efficiently. 

An approximate solution for the relic abundance is given by 

N(t)=NH = X ^ (3-12) 

. where x/=jr and Tf is the so called 'freeze-out' temperature and approximately, for 
the typical case x/ yy l,one has that Xf is a solution of 



y/xfe* f =y (3.13) 

We can then finally calculate the contribution of % to the energy density parameter 
finding a well know result 



Q y h 2 = 4x 1(T 10 (3.14) 

It is intriguing that so called "WIMPS"(e.g. the lightest supersymmetric particle) dark 
matter particles seem to reproduce naturally the right abundance since they have a 
weak cross section 



28 



0.1 



O.0L 



£ o.ooi 



LO- 



LO" 



lo: 



i I I I I r r r I i I F I I p r r I ■ 

; General Case of Comowia Number 

Density v J"* 


I ■ I I I r r i I I t I I'll 

\ Exponential drop 


Nonrelativistic Case of Comoving Number 
Density 


\\ Increasing < <™ > 


i i ■ i i i i i_i i i ■ i i ■ i i_i i 


\ -X : 

K .5 ; 

\ Freeze out: 



0.01 



0.1 



1 

x=in,T (time ■=£>-) 



10 



100 



Figure 3.1. 

Evolution of N and N eqn as function of x for m =100 GeV, g x = 2, g =90 (constant), 
o = 10~ 14 , 10" 1 Wl0~ 16 ,M Pi = 1.9x 10 l9 GeV 2 



a 2 tn aAOOGeV, 
a =^ = 10- 9 =— ) (3.15) 

and masses m x m ew lOOGeV. This observations is called the "WIMP miracle" and 
typically considered as an encouraging point supporting WIMPS as Dark Matter can- 
didates. 

The LW eqn. can be solved by approximately analytic calculation or numerical plotting 
by using Mathematica.We have plotted those equation by using both of the ways. As 
you can see Fig. 2.1. 
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LW Eqn., the termal production, is the most traditional mechanism but that many other 
mechanisms have been considered in the literature such as non thermal production of 
very massive particles(so called WIMPZILLAS) at preheating or even right-handed 
neutrino oscillations. 
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Chapter 4 



CONCLUSION 



The observational evidence for dark matter is overwhelming and now we are enter- 
ing one of the perhaps most exciting periods in physics history where a wide array 
of experiments may have the potential to unveil its nature. In this thesis I review some 
important concepts about dark matter. One of the best motivated candidates for dark 
matter so far is the WIMPS which I mentioned before. In this thesis the focus was 
on dwarf spheroidal galaxies, promising targets that are assumed to be strongly dark 
matter dominated systems due to their high mass-tolight ratios up to few hundred solar 
units. Dwarf galaxies of the Milky Way have the advantage of relative proximity and 
low astrophysical background as they do not lie within the galactic plane. Great impor- 
tance was given to studying the properties of these systems and how their dark matter 
profiles can be derived from kinematical studies. However a precise estimation of the 
density distribution is a challenging task and due to too little data available in particular 
the shape of the dark matter halo inner regions remains poorly known. Cusped profiles 
suggested from simulations like the Navarro-Frenk- White profile are not ruled out in 
general although observations seem to favor profiles with a flat inner core. 

The two investigated dwarf spheroidal galaxies Sagittarius and Canis Major are par- 
ticularly attractive targets as they are the closest dwarf galaxies discovered so far. The 
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uncertainty concerning their inner region is reflected in the plots of sensitivity limits 
for both investigated dwarf galaxies, where the distances between the two lines of a 
cored and a cusped profile are about one order of magnitude. 

Better opportunities for a significant dark matter detection could be obtained for in- 
stance by the next generation ground-based CTA (Cherenkov Telescope Array) exper- 
iment with bigger effective areas due to many Cherenkov telescopes in the order of 50 
and also IceCube experiment . Here, stereoscopy will be possible for small energies. 
Moreover the space based ray telescope FERMI will yield interesting results in the 
future. Apart from Sagittarius and Canis Major, there could be promising faint dwarf 
galaxies that are not discovered yet. However, there should be as many undiscovered 
dwarf galaxies in the other direction, too, which, when discovered, will provide inter- 
esting targets for the search of dark matter. 

In chapter 2 we gave a simple example of f(R) extended gravity in which the effect of 
DM is not due to exotic sources but due to the wrong choice of Lagrangian. Einstein's 
model, which extends that of Newton is characterized by the specific choice f(R) = R. 
In the f(R) model this has been modified to different polynomial powers of R. As an 
example the choice f(R) m R l i where «< 1, describes asymptotic behaviours of 
flat rotation curves and experimental results of Tully-Fisher (i.e. v 4 Mass). For n=0 
the model reduces to the Einstein's general relativity model which fails to explain flat 
rotation curves unless supplemented by exotic sources. 



32 



REFERENCES 



[1] M. Drees and H. Iminniyaz, Phys. Rev. D 76, 103524 (2007). 

[2] Kapteyn, Ap. J. 55, 302 (1992). 

[3] F. Zwicky, Helv. Phys. Acta 6, 110 (1933). 

[4] A. Slyz and T. Kranz, Mon. Not. Roy. Astron. Soc. 346, 1 162 (2003). 

[5] J. Bahcall and X. Fan, Astrophys. J. 504, 1 (1988). 

[6] R. G. Carlberg, Phys. Rept. 307, 67 (1998). 

[7] D. J. Eisenstein, Astrophys. J. 633, 560 (2005). 

[8] Hinshaw, Astrophys. J. Suppl. 148, 135 (2003). 

[9] NASA's Wilkinson Microwave Anisotropy Probe (WMAP), 



http://map.gsfc.nasa.gov/ 



[10] Particle Data Group, J. Phys. G. 33, 224 (2004). 

[11] G. Bertone, D. Hooper and J. Silk, Phys. Rep. 405, 279 (2005). 

[12] H. E. Haber and G. L. Kane, Phys. Rep. 117, 75 (1985). 

[13] WMAP Collab. and D. N. Spergel, Astrophys. J. Suppl. 170, 377, (2007). 

[14] E. W. Kolb, D. J. H. Chung and A. Riotto, [arXiv:hep-ph/98 1036T] 



[15] K. G. Begeman, A. H. Broeils and R. H. Sanders, MNRAS 249, 523 (1991). 

33 



[16] F. Zwicky, Ap. J. 86, 217 (1937). 

[17] R. H. Cyburt and H. Richard, Phys. Rev. D 70, 023505 (2004). 

[18] B. W. Lee and S. Weinberg, Phys. Rev. Lett. 39, 165 (1977). 

[19] A. S. Eddington, Cambridge University Press, Cambridge (1923). 

[20] H. Weyl, Ann. Phys. 59, 101 (1919). 

[21] S.Asgari and R.Saffari, Appl. Phys. Res. 2, 99 (2010). 

[22] T. P. Sotiriou and V. Faraoni, Rev. Mod. Phys. 82, 451-497 (2010). 

[23] S. Capozziello and M. De Laurentis, Annalen der Physik 524, 545 (2012). 

[24] S. Nojiri and S. D. Odintsov, arXiv: 0807.0685vl. 

[25] S. Nojiri and S. D. Odintsov, Phys. Rev. D 68, 123512 (2003). 

[26] S. Nojiri and S. D. Odintsov, Gen. Rel. Grav. 36, 1765 (2004). 

[27] R. Saffari and S. Rahvar, Phys. Rev. D 77, 104028 (2008). 

[28] P. Gondolo and G. Gelmini, Nucl. Phys. B 360, 145 (1991). 

[29] E. W. Kolb and M. S. Turner, Addison- Wesley (1990). 

[30] P. Di Bari, Some notes on Neutrino Physics and Particle Cosmology, 
http://152.78.192.50/^pdbld08/neupcimLpdf| ( 1 990) . 



34 



